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Introduction 


The National Aeronautics and Space Administration (NASA) is developing fast-time wake transport 
and decay models to safely enhance the capacity of the National Airspace System (NAS). These models 
are empirical algorithms used for real-time predictions of wake transport and decay based on aircraft 
parameters and ambient weather conditions. The aircraft dependent parameters include the initial 
vortex descent velocity ( V 0 ) and the vortex pair separation distance (b 0 ). The atmospheric initial 
conditions include vertical profiles of temperature or potential temperature (0), eddy dissipation rate 
(EDR), and crosswind. The model output consists of time history of vortex circulation strength and 
position. The fast-time wake models can be used for the systems level design of advanced air traffic 
management (ATM) concepts that safely increase the capacity of the NAS. It is also envisioned that at 
some later stage of maturity, these models could be used operationally, not only within the terminal 
airspace but also as onboard tools to support concepts such as dynamic separation of aircraft. 

NASA's first fast-time wake transport and decay model was developed by Greene (1986). In the late 
1990s, under NASA's Aircraft Vortex Spacing System (AVOSS) project, significant advances were made in 
wake vortex modeling based on the data from field experiments and large eddy simulations. The initial 
versions of the AVOSS Wake Vortex Prediction Algorithm (APA) were developed during the AVOSS 
program and demonstrated in the AVOSS demo at the Dallas/Ft. Worth (Hinton 2001). The APA model 
computes the out-of-ground-effect (OGE) decay and descent based on Sarpkaya (Sarpkaya 2000; 
Sarpkaya et al. 2001). The model has an algorithm for enhanced rate of decay during the ground effect 
developed by Proctor et al. (2000). The scheme to compute lateral vortex transport is based on the 
vertical profile of crosswind (Robins and Delisi 2002), and the in-ground-effect (IGE) transport accounts 
for vortex spreading and rebound (Robins et al. 2002). The code development of APA is described in 
Robins and Delisi (2002). 

NASA has also developed the TASS (Terminal Area Simulation System) Derived Algorithms for Wake 
Prediction (TDAWP) model. In the TDAWP model, the Sarpkaya component is replaced with algorithms 
developed from parametric studies using a Large Eddy Simulation (LES) model. The TDAWP model is 
described in Proctor et al. (2006) and Proctor and Hamilton (2009). The current version of the TDAWP 
model includes the effect of the crosswind shear gradient on transport. The TDAWP model is also 
denoted by TDP. 

The current distribution includes the latest versions of the APA (3.4) and the TDP (2.1) models. This 
User's Guide provides detailed information on the model inputs, file formats, and the model output. An 
example of a model run is also provided. A brief description of the Memphis 1995 Wake Vortex Dataset 
is also provided. Additional questions regarding this software distribution can be directed to Nashat 
Ahmad ( nashat. n.ahmad(5>nasa. gov) or Randy VanValkenburg ( randal.I.vanvalkenburggpnasa.gov ). 
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Software Distribution 

NASA's AVOSS Fast-Time Wake Prediction software distribution is given below: 

avoss/ 


' --bin/ 

I I 

| +-- apa34.exe 

I I 

| ' - - tdp2 1 . exe 


' --doc/ 

| '-- APA-UsersGuide.pdf 


' --etc/ 

| '-- cases . i 


' - -mem/ 

I I 

| +-- ADATA/ 

I I 

| +-- QDATA/ 

I I 

| +-- TDATA/ 

I I 

| +-- UDATA/ 

I I 

| +-- UPROXY/ 

I I 

| +-- CWP/ 

I I 

| '-- cws/ 

' --run/ 


Memphis 1995 Wake Vortex Dataset 

aircraft information and initial vortex location 
vertical profiles of eddy dissipation rate 
vertical profiles of potential temperature 
vertical profiles of crosswinds 
vertical profiles of proxy crosswinds 
lidar data for port vortex 
lidar data for starboard vortex 


+-- cases. i 

I 

+-- 1995-08-10-230029. apa34 

I 

+-- 1995-08-10-230029. tdp21 

I 

+-- 1995-08-10-230029. tplt 

I 

+-- 1995-08-10-230029. qplt 

I 

1995-08-10-230029. uplt 


list of cases to run 
APA3 . 4 output 
TDP2 . 1 output 

Potential Temperature (9) input 
Eddy Dissipation Rate (e) input 
Crosswinds input 
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Model Input/Output Filename Convention 

The file names are chosen to give a unique identifier for each case and model run. The unique 
identifier has the form: 


YYYY MO DY HRMNSC 


where, 



YYYY 

= 

Four digit year 

MO 

= 

Two digit month 

DY 

= 

Two digit day 

HRMNSC 

= 

Six digit HourMinuteSeconds 


Example: 

200 3_0 9_1 9_1 81019 has the following associated input files: 

200 3_0 9_1 9_1 81019. ADATA Initial Vortex Location & Aircraft Parameters 

200 3_0 9_1 9_1 81019. TDATA Vertical Profile of Temperature or Theta 

2003_09_19_181019 .UDATA Vertical Profile of Crosswind 

200 3_0 9_1 9_1 81019. QDATA Vertical Profile of Eddy Dissipation Rate 


The model output is written to (depending on the model used): 
2003_09_19_181019 . apa34 APA3.4 output 

200 3_0 9_1 9_1 81019. tdp21 TDP2.1 output 
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Model Input File Formats 

Case Input File (cases. i) 

The first six lines list the paths to input files. Line number 7 has the total number of cases (maximum 
number of cases = 5000) in the file and the type of lidar. The lidar type is not used in fast-time models - 
it is used by plotting and validation routines. The rest of the file lists the unique identifiers for each case. 
The path names should have a maximum length of 132 characters. 


CASE FILE (cases. i) 

/home / nnahmad/ apa -mode Is/ data / MEM1 995/ ADATA/ 

/home / nnahmad/ apa -mode Is/ data / MEM1 995/ QDATA/ 

/home / nnahmad/ apa -mode Is/ da ta /MEM1 995/ TDATA/ 

/home / nnahmad/ apa -mode Is/ da ta /MEM1 995/ UDATA/ 

/home / nnahmad/ apa -mode Is/ data / MEM1 995/ CWP / 
/home/nnahmad/apa-models/data/MEM1995/CWS/ 

2 0 ! total number of cases to run/lidar type 

1995-08-06-230412 

1995-08-06-232159 


Aircraft Parameters Input File ( YYY Y_M 0_D Y_H R M N SC . AD ATA) 

The first 11 lines of this file contain the header information. After the header lines, the file contains 
the initial vortex location and the aircraft data: Initial lateral location of the vortex (y 0 ), initial height of 
the vortex (z 0 ), initial descent velocity (V 0 ) i and the separation distance ( b 0 ) of the vortex pair. The units 
are MKS. 


AIRCRAFT PARAMETERS FILE (YYYY_MO_DY_HRM NSC. ADATA) 

# Location: MEM, 18L_TANG 

# Run Number: 1026 

# A/C Type: AT43 

# Wing span* (m) : 24.6 

# Weight* (kg) : 13940 

# ACspeed* (m/s) : 63.4 

# Air Density* (kg/m3) : 1.2 

# This file was created on 22-Mar-2012 12:42:37 Pacific Time 

# File created by Matt Pruis, NWRA, matt@nwra.com 

# *If data not provided with original data set, then default values used. 

# Data File Format: yo (m) , zo (m) ,Vo (m/s) , bo (m) 

5.2895, 90.03, 0.76635, 19.321 
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Temperature/Theta Profile ( Y YY Y_M 0_D Y_H R M N SC .TD ATA) 

The first 11 lines in this file contain header information. The first line after the header information 
gives the total number of data points in the file. If negative, then the file contains potential 
temperatures, otherwise the data are temperatures. The rest of the file has two columns: the first 
column lists the AGL height in meters and the second column lists the potential temperature in K or 
temperature in °C. If the user inputs a temperature profile, then it is converted by the model to 
potential temperature (9). 

At least three data points are required in the initial profile and the points should extend above and 
below the heights of vortex descent trajectory. If the input potential temperature profile contains 
regions of unstable stratification, then those values are set to zero (neutral stratification) within both 
the APA and TDP codes. Observations from various field sensors as well as simulation data from 
mesoscale models (Ahmad et al. 2013) can be used to generate the vertical temperature/potential 
temperature profiles. It should be noted that the fast-time wake models use potential temperature in 
model calculations. 


TEMPERATURE PROFILE DATA (YYYY_MO_DY_HRMNSC.TDATA) 

# Location: MEM, 18L_TANG 

# Run Number: 1026 

# A/C Type: AT43 

# Potential modifications from original data include : 

# (1) extrapolation above and below profile (with N=0 in these regions) 

# (2) removal of unstable regions that are not attached to ground 

# (with N=0 in these regions) 

# This file was created on 22-Mar-2012 12:42:37 Pacific Time 

# File created by Matt Pruis, NWRA, matt@nwra.com 

# Data File Format (1st line) : number of points (-ve indicates theta) 

# Data File Format (remainder of file) : z (m) potential temperature (K) 
-120 

0, 303.98 
5, 303.98 

10, 304.04 and so on 


5 



Crosswind Profile ( Y YY Y_M 0_D Y_H R M N SC . U D ATA) 

Similar in format to the xxx.TDATA file, but contains crosswinds (m/s). Heights are in meters. 
Crosswind profiles can be generated from the lidar data or estimated from the wake vortex trajectory 
(Pruis et al. 2011). 


CROSSWIND PROFILE DATA ( YYYY_M 0_D Y_H RM NSC. U DAT A) 

# Location: MEM, 18L_TANG 

# Run Number: 1026 

# A/C Type: AT43 

# Potential modifications from original data include : 

# (1) extrapolation above and below profile (with U constant in these regions) 

# This file was created on 22-Mar-2012 12:42:37 Pacific Time 

# File created by Matt Pruis, NWRA, matt@nwra.com 

# Data File Format (1st line) : number of data points 

# Data File Format (remainder of file) : z (m) Crosswind (m/s) 

120 

0, 2.124 

10, 2.074 and so on 


Eddy Dissipation Rate Profile (YYYY_MO_DY_HRMNSC.QDATA) 

Similar in format to the xxx . TDATA file, but contains the eddy dissipation rates ( m * 2 /s 3 * 5 * ). Heights are 
in meters. Given two observations of EDR, the vertical EDR profile can be generated using atmospheric 
boundary layer similarity theory (Han et al. 2000). EDR profiles can also be estimated from lidar data 
(Pruis et al. 2013). Input EDR values less than 10' 7 m 2 /s 3 are set to 10' 7 m 2 /5 3 within the models. 


EDDY DISSIPATION RATE DATA (YYYY_MO_DY_HRMNSC.QDATA) 

# Location: MEM, 18L_TANG 

# Run Number: 1026 

# A/C Type: AT43 

# Potential modifications from original data include : 

# (1) extrapolation above and below profile (with EDR constant in these regions) 

# This file was created on 22-Mar-2012 12:42:37 Pacific Time 

# File created by Matt Pruis, NWRA, matt@nwra.com 

# Data File Format (1st line) : number of data points 

# Data File Format (remainder of file) : z (m) EDR (m2/s3) 

120 

0, 0.0026156 

5, 0.0026156 

15, 0.002405 
20, 0.0023001 

30, 0.0020895 and so on 
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Lidar Data File Format 


Lidar Data File (YYYY_MO_DY_HRMNSC.cwp) 

For each case there are two lidar data files: one each for the port and the starboard vortices. The 
filename conventions are as follows: 


2 0 0 3_0 9_1 9_1 8 1 0 1 9 . cwp 
2 0 0 3_0 9_1 9_1 8 1 0 1 9 . cws 
200 3_0 9_1 9_1 8 1 0 1 9 . pip 
2003_09_19_181019.pls 


Continuous Wave Lidar Port Vortex Data 
Continuous Wave Lidar Starboard Vortex Data 
Pulsed Lidar Port Vortex Data 
Pulsed Lidar Starboard Vortex Data 


In most of the field experiments only one type of lidar was used and therefore a particular dataset 
will have either CW or PL files. In the Denver 2003 Field Experiment (Dougherty et al. 2004) both CW 
and PL lidars were deployed. 


The first six lines in the lidar data files contain the header information. Missing data points in the 
files are marked by a -9999. Please note that it is not uncommon for the lidar files to contain points 
that include valid times and vortex positions, but for which the circulation cannot be calculated due to 
missing or invalid data values. Consequently for a given vortex, the plots of circulation may display 
fewer discrete points than do the corresponding plots of lateral transport or altitude. 


LIDAR DATA FILE (YYYY_MO_DY_HRMNSC.cwp) 


# Location: MEM, 18L_TANG 

# Run Number: 1026 

# A/C Type: AT43 

# This file was created on 22-Mar-2012 12:42:38 Pacific Time 

# File created by Matt Pruis, NWRA, matt@nwra.com 

# Data File Format: time_p(s), y_pos_p (m) , z_pos_p (m) , Circ_p(m2/s) 
21 

8.88, 5, 87, 125.9091 
11.4, 8.3, 84.7, 158.1455 
13.6, 18.4, 92.5, -9999 

15.57, 20.5, 83.2, 123.7773 and so on 
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Model Output File Format 

Model Output File Format (e.g., YYYY_MO_DY_HRMNSC.apa34) 

The model output files are currently written in Tecplot® format. The model output file gives the 
time history of position (lateral distance and altitude) and circulation strength of both the port and 
starboard vortices. The output variables are listed in Table 1. 


Table 1: Fast-Time Wake Model Output 


column # 

variable 

description 

l 

time 

Time (seconds) 

2 

yp 

Lateral position of port vortex (m) 

3 

zp 

Altitude of port vortex (m) 

4 

gp 

Circulation of port vortex (m 2 /s) 

5 

ys 

Lateral position of starboard vortex (m) 

6 

zs 

Altitude of starboard vortex (m) 

7 

gs 

Circulation of starboard vortex (m 2 /s) 


The output filename extension is based on the model type: 

200 3_0 9_1 9_1 81019. apa34 APA3.4 output 

2003_09_19_181019 . tdp21 TDP2.1 output 

In addition to the model output, for each run the environmental initial conditions are also written to 
files for plotting purposes: 

2003_09_19_181019 .uplt Crosswinds 

2003_09_19_181019.qplt EDR 

2003_09_19_181019 . tplt Theta/Temperature 
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Model Output File Example (e.g., YYYY_MO_DY_HRMNSC.apa34) 

An example of a model output file is shown below: 


MODEL OUTPUT (YYYY_MO_DY_HRMNSC.apa34) 


TITLE="APA 3.4" 


VARIABLES = 

"Time (s) " 

, "Yp(m)", 

"Zp (m)", 

"Gp< 

;m A 2/ s) 

II 1 

r 

'Ys (m) " 

, "Zs (m) 

ZONE 

T="31" , 

. 1= 

= 1097 











0. 

.000 

8. 

.356 

41 

.821 

63. 

.793 

25. 

164 

41. 

.821 

63 

.793 

0. 

.100 

8. 

.300 

41 

.761 

63. 

.738 

25. 

108 

41. 

.761 

63 

.738 

0. 

.200 

8. 

.243 

41 

.700 

63. 

. 683 

25. 

051 

41. 

.700 

63 

. 683 

0. 

.300 

8. 

.187 

41 

. 640 

63. 

. 627 

24. 

995 

41. 

. 640 

63 

. 627 

0. 

.400 

8. 

.130 

41 

.580 

63. 

.572 

24. 

938 

41. 

.580 

63 

.572 

0. 

.500 

8. 

.074 

41 

.520 

63. 

.517 

24. 

882 

41. 

.520 

63 

.517 

0. 

. 600 

8. 

.017 

41 

.460 

63. 

.461 

24. 

825 

41. 

.460 

63 

.461 

0. 

.700 

7. 

. 960 

41 

.399 

63. 

.406 

24. 

768 

41. 

.399 

63 

.406 

0. 

.800 

7. 

. 904 

41 

.339 

63. 

.351 

24. 

712 

41. 

.339 

63 

.351 

0. 

. 900 

7. 

.847 

41 

.279 

63. 

.296 

24. 

655 

41. 

.279 

63 

.296 

1. 

.000 

7. 

.790 

41 

.220 

63. 

.241 

24. 

598 

41. 

.220 

63 

.241 

1. 

.100 

7. 

.734 

41 

.160 

63. 

.186 

24. 

542 

41. 

.160 

63 

.186 

1. 

.200 

7. 

. 677 

41 

.100 

63. 

.131 

24. 

485 

41. 

.100 

63 

.131 

1. 

.300 

7. 

. 620 

41 

.040 

63. 

.076 

24. 

428 

41. 

.040 

63 

.076 

1. 

.400 

7. 

.563 

40 

. 980 

63. 

.021 

24. 

371 

40. 

. 980 

63 

.021 

1. 

.500 

7. 

.507 

40 

. 921 

62. 

. 966 

24. 

315 

40. 

. 921 

62 

. 966 

1. 

. 600 

7. 

.450 

40 

.861 

62. 

. 912 

24. 

258 

40. 

.861 

62 

. 912 

1. 

.700 

7. 

.393 

40 

.802 

62. 

.857 

24. 

201 

40. 

.802 

62 

.857 

1. 

.800 

7. 

.336 

40 

.742 

62. 

.802 

24. 

144 

40. 

.742 

62 

.802 

1. 

. 900 

7. 

.279 

40 

. 683 

62. 

.748 

24. 

087 

40. 

. 683 

62 

.748 

2. 

.000 

7. 

.222 

40 

. 623 

62. 

. 693 

24. 

030 

40. 

. 623 

62 

. 693 


and so on 
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Memphis Field Experiment (1995) 

A comprehensive field experiment to measure wake vortices and the associated ambient 
meteorological conditions was conducted at the Memphis International Airport in Memphis, Tennessee 
from August 6 through August 29, 1995 (Zak 1995; Campbell, et al. 1997). The experiment was 
sponsored under NASA Langley Research Center's Aircraft Vortex Spacing System (AVOSS) project 
(Hinton 1995; Perry et al. 1997). The wake data were collected using a continuous wave lidar (Figure 1). 
The meteorological sensors included radiosondes, sodars, a wind profiler, one 150 ft high meteorological 
tower, a Radio Acoustic Sounding System (RASS), and NASA Langley's OV-10 research aircraft. The 
radiosondes were used to measure winds and temperature measurements (10s averages) at 50m 
vertical resolution. The OV-10 aircraft was flown at selected times and took measurements of 
temperature and winds at a sample rate of 10 Hz. Temperature (5 min averages) was measured using 
RASS every 30 min at 14 vertical levels from 127m to 1492m. The 150/f (45.7m) meteorological tower 
was equipped with a large array of sensor systems. Winds, temperature and moisture were measured 
from the tower at 5m, 10m, 20m, 30m, and 42m heights. Turbulence quantities (turbulence kinetic 
energy and eddy dissipation rate) were estimated from wind measurements at 5m and 40m heights. 
Rain rate, soil temperature, soil moisture, barometric pressure, and incoming and outgoing solar 
radiation also were measured by the sensors deployed on the meteorological tower. Standard 
meteorological data such as atmospheric pressure, temperature, moisture, cloud cover, visibility, etc. 
were obtained from the National Weather Service's Surface Aerodrome Observations (SAO) and the 
Automated Surface Observations System (ASOS). 

Data Processed for Fast-Time Wake Models 

Aircraft Data 

The aircraft data used by the fast-time wake prediction models include the initial position (offset) of 
the vortex pair with respect to the runway centerline (y 0 ), the initial height of the vortices (z 0 ), the initial 
vortex descent rate ( V 0 ) and the initial separation of vortices (b 0 ). 

The initial position (offset) of the vortex pair with respect to the runway centerline y 0 was estimated 
using an average of the first few data points for each landing. The initial height of the vortices z 0 was 
estimated from backward extrapolation of the altitude trajectory in time. The initial separation distance 
between the vortices b 0 was estimated assuming an elliptical wing loading, 



(1) 


where B is the wingspan of the aircraft. The initial vortex descent rate was estimated from the aircraft 
weight, aircraft speed, air density, and the initial vortex separation b 0 , 


gfV 

InpUbl 


( 2 ) 


where g (=9.81m/s 2 ) is the acceleration due to gravity, p is the air density - which was assumed to be 
1.2kg/m 3 for all the landings, U is the reported airspeed, and W is the reported landing weight of the 
aircraft. 

Types of aircraft observed at different measurement sites are listed in Table 2 and a graphical 
presentation of aircraft distribution is shown in Figures 1-2. The Armory site was located south of the 
airport, while the TANG, Tchulahoma, and Threshold sites were all located at the north end of the 
airport. 
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Crosswinds and Headwinds 

The profiles of the mean crosswind and headwind were generated by the AVOSS Winds Analysis 
System (AWAS) using an optimal estimation of data fusion from several different wind sensors including 
two Doppler radars, the meteorological tower, and the SODAR. There are some known deficiencies in 
the AWAS profiles which are discussed in detail by Dasey et al. (1998). In the current distribution, only 
the crosswinds are included. 

Eddy Dissipation Rate 

EDR profiles were estimated using the two sonic anemometers on the meteorological tower and 
extrapolating to heights using atmospheric boundary layer similarity theory (Han et al. 2000). The 
profiles were generated using the atmospheric turbulence profile generator (ATPG) code which 
implements the algorithms described in Han et al. (2000). Turbulence profiles were extrapolated to the 
ground (z = 0) and to a height above the observed vortices with a constant EDR value whenever the 
measured profiles did not extend to those heights. 

Stratification 

Temperature profiles were estimated using a fusion of the RASS and temperature sensors on the 
ASOS and the meteorological tower. The temperature profiles were converted to potential 
temperatures using the dry adiabatic lapse rate. The highest temperature measurement on the tower 
was at approximately 43m AGL and the lowest observation of the RASS was at 127 m. A known 
deficiency in the profiles is a frequent mismatch in these two temperatures, leading to unrealistic 
gradients in the temperature profile within this region. This can sometimes lead to highly unstable 
persistent regions that are not attached to the ground. The potential temperature profiles were 
therefore pre-processed to remove these unstable regions by making the potential temperature 
constant in these regions. 


Table 2: Aircraft Observed at different Measurement Sites 


Aircraft Type 

Armory 

TANG 

Tchulahoma 

Threshold 

Total 

AT 4 2 

1 

1 

- 

- 

2 

B727 

86 

1 

3 

24 

114 

B7 37 

3 

- 

- 

- 

3 

B7 57 

5 

2 

1 

- 

8 

BA31 

- 

- 

2 

- 

2 

DC10 

25 

- 

- 

2 

27 

DC 9 

57 

19 

6 

3 

85 

EA30 

11 

- 

- 

- 

11 

EA31 

9 

- 

- 

1 

10 

EA32 

15 

7 

2 

1 

25 

FK10 

5 

2 

- 

- 

7 

MD1 1 

- 

- 

- 

1 

1 

SF34 

2 

- 

8 

- 

10 

Total 

219 

32 

22 

32 

305 
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Figure 1: Distribution of the aircraft types in the MEM dataset is shown in the left panel. A B757 
landing in the background of the NASA lidar van is shown in the right panel. 


250 m 





4000 m 

Figure 2: Distribution of the lidar data by aircraft category and location. The Armory site was located 
south of the airport, while the TANG, Tchulahoma, and Threshold sites were all located at the north 
end of the airport. 
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Evaluation of Fast-Time Models using the Memphis 1995 Data 

Several evaluations of NASA's fast-time models have been conducted in the past (Proctor and 
Hamilton 2009; Pruis and Delisi 2011a; Feigh et al. 2012). Under a NASA Research Announcement 
(Enabling Super-Dense Operations by Advancing the State of the Art of Fast-Time Wake Vortex 
Modeling), the Northwest Research Associates (NWRA) was tasked to conduct an independent 
evaluation of NASA's fast-time models over a three year period. This evaluation concluded that in 
general the errors in model circulation predictions had a mean root mean square error on the order of 
0.2r 0 to 0.3r 0 (r 0 is the initial wake circulation), the vertical transport errors were on the order of 0.5b 0 
and the lateral transport errors were on the order of b 0 (Pruis and Delisi 2011b). NWRA also 
demonstrated that the lateral transport errors can be reduced to as low as 0.5b 0 if more accurate 
crosswind initial conditions (e.g., by using proxy crosswinds as initial conditions) were provided to the 
fast-time models (Pruis et al. 2011). 

In this section the current distribution of the fast-time models are evaluated using the continuous 
wave lidar observations from the Memphis 1995 field experiment. The accuracy of predictions for the 
two models was quantified in terms of root mean square error (Error rms ), mean absolute error {Error mae ), 
and Bias: 


Error 

rms 



Error 

mae 



-x° bs \; Bias = — V 
' n 



(3) 


The prediction errors in TDP2.1, and APA3.4 for all Memphis cases are given in Table 3. The errors in 
TDP2.1 and APA3.4 categorized approximately by phase (OGE, NGE, and IGE) are given in Tables 4-6. 
Proxy crosswinds were not used in this evaluation. 

For the Memphis data, the groupings of the landings were based on the measurement site where the 
observations were obtained. The OGE observations are the measurements obtained at the 36R_Armory 
site, where the mean initial vortex height was 180m. The NGE observations were obtained from the 
18L_TANG site where the mean initial observations, were at a height of 99m. The IGE data came from 
two different lidar locations; the 27_Tchulahoma and the 27_Threshold data. The mean height of the 
initial vortex observation for these two sites was 36m. In dimensional space, the combined different 
phases of vortex trajectory (OGE, NGE and IGE) are defined as follows: 

z o >130»; => OGE 

75»; < z a < 130m =4> NGE (4) 

z 0 < 15m => IGE 


Table 3: Memphis: All 305 Cases 


Model 

Circulation 
(normalized by r 0 ) 

Lateral Transport 
(normalized by b 0 ) 

Altitude 

(normalized by b 0 ) 


rmse 

mae 

bias 

rmse 

mae 

bias 

rmse 

mae 

bias 

TDP2.1 

0.26 

0.22 

0.02 

1.01 

0.83 

0.09 

0.52 

0.44 

0.07 

APA3.4 

0.24 

0.21 

- 0.05 

0.97 

0.80 

0.10 

0.54 

0.46 

0.14 


13 



Table 4: Memphis OGE: 219 Cases 


Model 

Circulation 
(normalized by r 0 ) 

Lateral Transport 
(normalized by b 0 ) 

Altitude 

(normalized by b 0 ) 


rmse 

mae 

bias 

rmse 

mae 

bias 

rmse 

mae 

bias 

TDP2.1 

0.25 

0.21 

0.03 

1.06 

0.87 

0.16 

0.58 

0.50 

0.06 

APA3.4 

0.23 

0.20 

- 0.04 

1.02 

0.83 

0.18 

0.60 

0.51 

0.15 


Table 5: Memphis (NGE): 32 Cases 


Model 

Circulation 
(normalized by To) 

Lateral Transport 
(normalized by b 0 ) 

Altitude 

(normalized by b 0 ) 

rmse 

mae 

bias 

rmse 

mae 

bias 

rmse 

mae 

bias 

TDP2.1 

0.25 

0.21 

- 0.02 

1.05 

0.86 

- 0.32 

0.46 

0.39 

0.12 

APA3.4 

0.25 

0.21 

- 0.12 

1.04 

0.86 

- 0.33 

0.47 

0.39 

0.18 


Table 6: Memphis (IGE): 54 Cases 


Model 

Circulation 
(normalized by T 0 ) 

Lateral Transport 
(normalized by b 0 ) 

Altitude 

(normalized by b 0 ) 

rmse 

mae 

bias 

rmse 

mae 

bias 

rmse 

mae 

bias 

TDP2.1 

0.31 

0.26 

0.02 

0.67 

0.58 

0.02 

0.23 

0.20 

0.06 

APA3.4 

0.31 

0.26 

- 0.05 

0.67 

0.57 

0.01 

0.22 

0.20 

0.02 
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Example from the Memphis 1995 Dataset 

An example case is given in the run/ directory of the distribution. The directory, contains the 
following cases . i file: 


cases, i 

/ home / nnahmad/ apa -mode Is/ avo s s /mem/ADATA/ 
/home /nnahmad/apa -model s/avoss/mem/QDATA/ 

/ home / nnahmad/ apa -mode Is/ avo s s /mem/ TDATA/ 
/home /nnahmad/ apa -model s/avoss/mem/UPROXY/ 
/home /nnahmad/ apa -model s/avoss /mem/ CWP/ 
/home /nnahmad/ apa -model s/avoss/mem/CWS/ 

1 0 

1995 - 08 - 10-230029 


The user needs to edit the first six lines to reflect the location of the input files. In this case, the 
proxy crosswinds were used. The number of cases is 1 and the lidar flag is set to 0 for the CW lidar. The 
lidar flag is not used by the fast-time models and is used only by post-processing routines. After editing 
the cases, i file, the fast-time models tdp21.exe and, apa34.exe can be invoked. Figure 3 shows 
the environmental initial conditions as well as the * .ADATAfile for this case. The output of APA3.4 and 
TDP2.1 is compared with the CW lidar data in Figures 4-5. 


Environmental Initial Conditions 

ADATA File 



theta (K) 



300 302 304 306 308 310 



-5 

u (m/s) 
0 

5 

# Location: MEM, 36R Armory 

600 

1 1 

; Case ID 

1995-08-10-230029 

1 1 1 


# Run Number: 1108 




/ 



# A/C Type: DC93 

500 



/ 



# Wing span* (m) : 28.5 







# Weight* (kg) : 46653 

400 






# ACspeed* (m/s) : 78.1 

e 






# Air Density* (kg/m3) : 1.2 

N 300 






# File was created on 22-Mar-2012 







# File created by Matt Pruis , NWRA 

200 













# 

100 






# yo (m) , zo (m) ,Vo (m/s) , bo (m) 







17.301, 218.589, 1.5512, 22.384 

fc 


..ml . . 




>-* 10‘ 7 10"* 

10' ! 

10"‘ 10- 3 




edr (m A 2/s A 3) 




Figure 3: The environmental initial conditions (crosswind, eddy dissipation rate, and potential 

temperature are plotted in the left panel. The ADATA file for this example case is given in the right 
panel. 
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Figure 4: The circulation strength ( m 2 /s ) predicted by APA3.4 and TDP2.1 are plotted along with the 
lidar data in the left panel. The right panel shows the comparison of vortex descent with lidar 
observations. 




Figure 5: The vortex trajectories predicted by APA3.4 and TDP2.1 are plotted along with the lidar 
data. 
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